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ABSTRACT 

A  self-contained  treatment  of  the  theory  of  radio  waves  in  an 
inhomogeneous  atmosphere  is  given.     The  refractive  index  is  assumed 
to  vary  with  height  above  the  earth  in  a  monotonic  fashion.     Variation 
according  to  an  exponential  law  is  used  for  illustration  of  general  prin- 
ciples.    For  this  case,    rigorous  series  formulas  are  developed  for  the 
distance  to  the  horizon  for  an  elevated  point  in  the  atmosphere. 


ON  RADIO  WAVE  PROPAGATION  IN  AN 
INHOMOGENEOUS  ATMOSPHERE 


1.     Introduction  and  Formulation 

It  is  the  purpose  of  this  paper  to  review,    and  extend,    work  of 

1  2  3  4 

Watson    ,   Eckersley  and  Millington    ,   Bremmer    ,    Vvedensky     and 

5 
Fock     on  propagation  of  radio  waves  in  an  inhomogeneous  atmosphere. 

In  particular  an  attempt  will  be  made  to  present  a  unified  treatment  which 

includes  the  analyses  of  the  above  authors  and  shows  that  'Ithe  distance 

to  the  horizon"  is  a  valid  concept  in  diffraction  theory  even  when  the 

atmosphere  is  not  homogeneous.     The  earth  is  to  be  represented  by  a 

homogeneous  sphere  of  radius  a  and  the  source  is  to  be  a  vertical 

electric  dipole  which  for  the  moment,    is  located  near  the  earth's 

surface.     The  usual  spherical  coordinates  (r,    9,   0)  are  chosen  with 

the  source  oriented  along  the  polar  axis  and  the  surface  of  the  earth 

is  defined  by  r  =  a.     The  dielectric  constant  relative  to  free  space  is 

denoted  by  K  and  is  understood  to  be  a  function  of  r  and  decreases 

monotonically  with  increasing  r. 

The  resultant  magnetic  field  has  only  a  single  component,   H  ,, 

because  of  the  assumed  polar  symmetry  in  the  problem.     As  a  further 

consequence  of  the  symmetry,   the  electric  field  component  E  .  is  zero. 

It  is  convenient  to  express  the  remaining  field  components  in  terms 

of  single  scalar  function  U  which  is  actually  the  radial  component  of 

a  Hertz  vector.     For  a  time  factor  exp  (iojt),   this  leads  to 

H0  =  ik°  K  -Vi  '  <*> 


_  120rr  9      /   .     a      9U\  , 

E       =    = — 7, —     -K-Tp  (  sin  9     -5-75-    )  ,    and 

r  r  sin  9  9  9    V  d  9      ' 


(2) 
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E9    =   XT     Tr.  (K  r   1T0  <3> 

where  k0     =     Zit/wavelength  in  vacuo. 


-   2  - 


These  field  components  satisfy  Maxwell's  equations  if  the  function  U 
satisfies 


uV[4-  ^(Kru» 


r  sinfl    ~d!T 


sin0    -|y-    1    -J   k!  E    i   U     =    © 


Following  the  work  of  Foe k    ,    it  now  proves  to  be  convenient  to 
introduce  an  auxiliary  quantity  U1  which  is  defined  by 

U»     =    e"lka0    K  r     N/sin  0      U 
and  consequently  the  equation  for  U  satisfies 


(4) 
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where  K'     =    dK/dh,   K     is  the  value  of  K  at  r  =  a,   and  k  =  k.    n/k" 

o  o 

is  the  effective  wave  number  at  the  surface  of  the  earth. 


Since  E 


0 


-  ZH,  at  r  =  a    it  follows  that 
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au1 
a  r 


ik  (K   )2  ZU'     at  r  =  a 
o 


(7) 


6-8 


where  Z  is  the  effective  value  of  the  surface  impedance 

The  right-hand  side  of  equation  (6)  can  be  replaced  by  zero.     The 
justification  for  doing  this  is  based  on  the  fact  that 


9U'/9  r      =    Order  of  (kU'/M) 


and 


9U'/90    =    Order  of  (kaU'/M   ) 

1/3 
where  M  =  (ka/2)   '        is  a  large  quantity.     It  is  also  possible  to 

simplify  equation  (  6  )  by  noting  that 


(1 


!/r2)     =*     2h/a    . 


-  3  - 
Consequently  U'  satisfies 

where  s  =  ad  is  the  great-circle  distance  measured  along  the  earth's 
surface.     The  boundary  condition  at  the  surface  of  the  earth  then 

becomes 

l 

8U'/9k    =    +ik(K)2ZU!    (at  h  =  0)      .  (9) 

The  function  U',    which  must  satisfy  equations  (8)  and  (9)» 

should  also  have  the  proper  singularity  at  the  source.     That  is  it  should 

behave  as 

-ikR 

U    =   ^~    (1  +Rv)  (10) 

near  the  source,    where  R     is  a  Fresnel  reflection  coefficient  of  the 

v 

ground  for  vertically  polarized  waves.     For  h  <  <  s,    and  assuming 
that  our  approximate  boundary  condition  is  valid,   R  s  s  +  (h2/2s) 

and  thus  Rg    ffi;g    .  (11) 

Consequently,   the  required  behavior  of  U'  near  the  origin  is 

kh2 
i  -1  -i  -t — 

U'^K     (a/s)2    (2h/s)  (h/s  +  Z)    e         S  (12) 


2.     Equivalent  Earth  Radius  Concept 
The  equivalent  earth  radius  a     will  now  be  defined  in  the  following 


manne  r 
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where 

K' 
o 

= 

™\  =  o   ■  (13> 

This  equation  is  usually  justified  on  the  basis  of  geometrical  optics. 
It  can  also  be  demonstrated  in  a  more  rigorous  fashion  by  rewriting 
equation  (8)  in  the  form 


l^-Zik-W+^lhd+tfU.    =0  (14) 
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Wh6re  ae        /  K  "  Ko 

g  =   -2K~    ( ET~^    "    Ko 


Now  if  the  atmosphere  is  linear  then  by  definition    K'     =    (K  -  K   )/h 
and  y  =  0  and  consequently 


82U'  -,„_  8U1      .      2k2 

a 


W~   "  2ikTTi     +  4^  hU'  =  0  (15) 


But,    on  the  other  hand,   if  K  =K     for  all  h  then  from  equation  (8)  or 
(14), 

B2  U»  -,„_   9  U1      .      2k2 

a 


STpr-   "    akTT-  +   ^~    h  U'   =   °     ■  (16» 


This  strongly  suggests  that  propagation  in  a  linear  atmosphere  of 
radius  a  is  equivalent  to  propagation  over  an  earth  of  modified 
radius  a     with  a  homogeneous  atmosphere.     It  is  therefore  concluded 
that  the  effective  earth  radius  is  indeed  a  valid  concept.     This  con- 
clusion has  also  been  reached  by  many  other  authors,   including 
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Bremmer,   Fock  and  Miller    . 

3.     Extension  to  Non-Linear  Atmosphere 
In  many  practical  cases  of  interest,  however,  K  is  not  a  linear 
function  of  height,  particularly  at  great  heights.     Then  the  function  g 
is  non-zero  and  is  itself,   in  general,    a  function  of  height.     To  further 
analyze  this  case,   the  revelant  equations  will  be  put  in  dimensionless 
form.     The  following  parameters  are  first  introduced 


i 


i 

1 
/  ka    \  3  1 

-i  (_e   \       (K  j  3    (z/120tt),    and 
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1/6 


w   =  IT"  (-TTI-  ■ 

o       \  < 
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The  equation  for  W  is 

_|^L     .      i|it»|Hg|W    =    0  (17) 

and  the  boundary  condition  at  the  earth's  surface  now  becomes 

-JPL    +    qW    =    0aty    =    0.  (18) 

Furthermore,    as  x  tends  to  zero,    W  should  behave  as  (2/\Tx)  • 
exp  (-iy2/4x). 

It  is  now  desirable  to  first  consider  the  solution  of  equation  (17) 
when  y  can  be  neglected  compared  to  unity.     This  corresponds  to  an 
atmosphere  which  is  linear  or  near  linear.      The  revelant  equation  is 
now 

a2  w  aw 


sw 


-   i-xir  +   y  w  =  °   •  (19) 


ix 


This  can  be  solved  by  the  method  of  separation  of  variables.     That  is, 
assuming 

W    =   X  (x)    Y(y) 
one  finds 

12 


1  d^Y(y)  _      .    _1  dX(x)    _ 


(20) 


where  t  the  separation  parameter.     Thus 

JET    +    itx    =    0  (21) 

and 

4^-    +    (y  -  t)  Y    =    0    .  (22) 

Solutions  of  the  first  equation  are 

X  (x)     =    e±itx  (23) 

and  the  second  equation  has  solutions  of  the  form 

Y(y)     =    W(t  -  y)  (24) 

where  W(t)   satisfies  Airy' s  equation 
d2W(t) 


"oT^ 


t  W(t)     =    0     .  (25) 
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The  appropriate  solution  for  the  x  dependence  is  chosen  such  that 
outgoing  waves  are  only  retained.     Therefore 

X(x)    =    e"ltx  .  (26) 

Furthermore,   the  radial  part  of  the  solution  (i.  e.    y  dependence)  is 
chosen  so  that  the  phase  of  Y(y)  is  a  decreasing  function  with  height 
(i.e.   phase  lag  increases  with  height).     At  this  stage,    it  should  be 
noted  that  Airy's  equation  is  closely  related  to  Bessel's  equation  of 
order  1/3.     In  fact 

I  ?A 

W(t)    =    const,  x  t2    Zl  /3  [(2/3)  (-t)   '    ]  (27) 

where  Z .  /_  is  any  cylindrical  Bessel  function  or  linear  combinations 

thereof.     In  order  to  satisfy  our  requirement  that  the  phase  lag  be  an 

(2) 

increasing  function  with  height,   the  Hankel  function  H.  /~  of  the  second 

must  be  used.     A  particular  solution  which  is  suitable  from  the  above 
standpoint  is 

W(t)    =    Wx(t)    =    e"i2ir/3    (-ut/3)2    H(^3[(2/3)(-t)3/2]  (28) 

where  the  constant  is  chosen  to  be  consistent  with  the  work  of  Fock  and 
Leontovich.     In  fact,   the  latter  authors  represent  W,(t)  as  an  integral 
of  the  Airy  type 

l  C°  t  3/3 

»1«    =     ^     Je12,/3      *       '  dZ  «29» 

where  the  contour  is  a  straight  line  from  infinity  along  the  straight 
arg.  t  =  2tr/3  to  the  origin  and  then  out  to  oo  along  the  real  axis, 
(i.e.   arg.   t     =    0). 

To  satisfy  the  boundary  condition  at  the  earth's  surface,    one  must 
choose 

[W'jft-y)     -    qWjIt-y)]  =    0     .  (30) 

A  numerical  study  of  this  equation  indicates  that  the  roots  t  =  t 
(s  =  0,    1,    2.    .    .)  are  in  the  fourth  quadrant  of  the  t  plane.      The 
solution  is  then  made  up  a  linear  combination  of  terms  of  the  type 

—  lXt c       ,.,»  , 

Wjft8-y). 
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It  can  be  seen  immediately  that  an  equivalent  representation  is 

r>       e_ixt    Wj(t-y) 

Wl     =      V    J  W'    (t)     -     q  W.(t)       dt  <31> 

C 

where  C  is  some  contour  which  encloses  the  poles  of  the  integrand 
(i.e.   at  t  =  t  ).     The  factor  y,   assumed  constant,   is  to  be  chosen 
such  that  W,  has  the  proper  singularity  as  x  approaches  zero.     It  is 
to  be  expected  that  behavior  of  W,  for  small  x  is  to  be  determined  by 

the  behavior  of  the  integrand  for  large  values  of  t  because  of  the 

-ixt 
factor  e         »     But,   when  t  is  very  large,   it  follows  from  the  known 

asymptotic  behavior  of  the  Hankel  functions,   that 

W,  (t-y)  1  i 

W[(t)  -  q  Wjt)      "      ±    (t)"2exp[-y(t)2]  (32) 

where  the  upper  (+)  sign  is  to  be  employed  for  2ir/3    =  arg  t  >  -tt/3 
and  the  lower  (-)  sign  for  -ir/3  >  arg  t  ^  4tt/3.     The  contour  C  is  now 
drawn  such  that  one  branch    is  along  the  negative  imaginary  axis 
(-  i  oo  to    0)  and  the  other  along  the  real  axis  (0  to  oo).     For  small 
values  of  x,   then 


C  -ixt     e^ 

W  ~  v  \  e  en_        dt  (33) 


-ysTT 
+    v     \       e"""'   -1 dt 


f         -ixt      e 


nTT 
These  integrals  can  be  transformed  to  the  form 


00 

j    e"p2  dp  =  *nr/z 


(34) 


and  thus  l^         .    y2 

w  ~  2y  (-s— )"  e     **  (35) 

1 
and  therefore,    if    y   =  (i/tr) 2 
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W   as    2  x"2    e  (36) 


as  x  approaches  zero. 

The  required  solution  can  thus  be  written 

e~ixt  W.(t-y) 
W 


=  (77     Jc   Wj(t)  -  q  w(t)     dt   •  <37> 


In  this  particular  case  the  solution  for  U  becomes 

-i* 

s 


-ika0  -ika0 

U    =   — i    V    s    -£-, V  (38) 

a(0  sin  0)2  ay 

where  V  =  x2  W 

/ixv  r    e" xt  wi<*-y) 

V    =    ("T  )       Jc      WF>  -  q  WlW    dt    •  <39> 

Attention  is  now  returned  to  the  case  of  an  inhomogeneous  atmosphere 
wherein  g(y)  is  no  longer  negligible  compared  to  unity.     The  more 
general  equation  for  W  must  then  be  retained 

4jr~  '    *  -%T  +   y  [  l  +  g(y)]  w  =  °  '  (40> 

Proceeding  in  an  analogous  fashion  as  in  the  linear  case,   it  is  readily 
found  that  solutions  are  of  the  form 

W    =    e"lxt    F  (y,t)  (41) 

where  F(y,  t)  must  satisfy 

d2^Y,t)       +    [y-t    +    yg(y)]   F(y,t)    =    0    .  (42) 


The  solution  is  then  expressed  as  an  integral 

•c    r-^-  +  qrn 
L  J  v=o 
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which  is  constructed  so  the  poles  of  the  integrand  are  the  roots, 


t  =  t    ,    of  the  equation 


9F(y,t) 
9y 


+     qF(y,t) 


=  0 


(44) 


y=0 


and  the  constant  term  outside  the  integral  is  chosen  so  that  it  reduces 
to  the  corresponding  expression  for  W,   in  the  linear  atmosphere 
(i.e.    g  =  0).     The  corresponding  (formal)  solution  for  the  function  U 
is  thus  written  .,     n 

U    = i      V 


where 


V    = 


H-li 


a  (0  sin  0)2 

-ixt  _.      ,. 
e  F(y,t) 

L.  -Jv=o 


dt 


(45) 


(46) 


4.     Asymptotic  Form  of  the  Solution 

Various  asymptotic  approximations  are  now  investigated.     For 

large  values  of  the  real  part  of  (y-t    +    yg),   the  function  F  (y,t)  of 

equation  (40)  is  given  approximately  by  the  W.K.B.    solution  of 

equation  (40).     This  is 

.-iir/4 


F(y,t)     a 


Ce 


[y-t    +    yg  (y)]4 


exp 


_-i  r  [  u- 

L     Jt 


t    +    ug(u)]2  du 


(47) 


where  C  is  a  constant  and  the  lower  limit  of  the  integral  is  arbitrary. 
It  should  be  noticed  that  if  g  =  0 

F(y,t)    =    C  Wx(t-y) 
which  for  large  t-y  behaves  asymptotically  like 


F(y,t) 


C  e 


-iir/4 


exp 


[-{^ 


du 


(48) 


(y-t)  4 
in  accordance  with  equation  (28). 

Equation  (46)  for  V  in  the  case  of  an  inhomogeneous  atmosphere 
can  now  be  simplified  by  noting  F(y,t),    where  it  occurs  in  the  denominator, 


can  be  replaced  by  W,  (t-  y)  and  therefore 


3F 


+    qF 


> 


y  =  0 


W'x(t)  +  qWj(t). 
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In  the  numerator  F  can  be  replaced  by  the  asymptotic  form  given  by 
equation  (46)  where  y  is  large.     Thus 


V    =    - 


where 


1    /  x 


T  II? 


(0i 


1 

y  -  t    +  yg(y) 


-ifl 


WJ(t)  -  qWjtt) 


dt        (49) 


Q    =    xt    +      \         n/u  -  t    +    ug(u)     du      . 
Jo 


The  saddle  point,    determined  by  d^/dt  =  0,    is  thus  given  by  t  =  -p2 
where  y 

f  du  -    2x    =    0  .  (5 


Jo       ^/u  +  p2     +    ug(u) 


Thus,   following  the  prescription  of  the  saddle  point  method 


-too 


2x 


8p 


where 


and 


F( 


"p)  =  3T  I 


ipt 


F  (-p) 


dt 


W»(t)   -  qWj(t) 


ftc    =    -xp' 


it- 


+  p2     +    ug(u)]  2    du 


(51) 
(52) 

(53) 


The  integral  F(-p)  has  been  discussed  on  previous  occasions  and 

10 
numerical  results  are  available      .     For  large  negative  values  of  the 

argument  (i.e.   p   >  0)  it  was  noted  that 


F(-p)  =    plP{q     exp.(-ip3/3) 


(54) 


and  therefore 


e-iJi.    f2x   _^P.^2  Zp;  exp.(_ip3/3) 


By     /  p  +  iq        ^"^"»    "r   '  "'  *55* 

must  be  the  geometrical  optical  approximation  and  consequently,   the 
factor  2p/(p  +  iq)  can  be  identified  as  1  +  R(a)  where  R(a)  is  the  Fresnel 
reflection  coefficient  for  a  plane  wave  incident  at  an  angle  a  on  the 
boundary.     The  parameter  p  is  then  related  to  a  by 


-  11  - 

1 

p=  (ka  /2)3  cos  a  . 

In  the  case  of  the  linear  atmosphere  (i.  e.    g(u)  =  0) 

fl.    =    -xp2    +    \      (u    +    p2)  2  du  (56) 

Jo 

m  p3/3    +    2y3/2/3 
for  p  small  compared  to  x  and  y.     In  this  same  case 


2 


P    =  -^^  <57) 


and  consequently  i_ 

2 

£*     1 


(-£") 


When  g(u)     =    0,   the  formulas  thus  reduce  to  those  for  a  homogeneous 
atmosphere  with  a  modified  earth  radius. 

In  the  case  of  the  inhomogeneous  atmosphere  equation  (50)  for  p 
can  be  solved  approximately  to  yield 

p    =  4-     \  d"  -  x  (58) 

C  ^0  N/    U     +      Ug(u) 

under  the  condition  that  x  and  y  are  large  but  the  difference  x  -  NTy" 

is  finite.     Thus 

x  \ 

(**  ^r)  '      =    ~ r        =     1  (59) 

V      dy  )  [y  +  yg(y)l* 

and  y 

Q0     &  p3/3    +      V     [u  +  ug(u)]2    du  .  (60) 

Jo 

For  application  to  an  actual  atmosphere,   it  is  necessary  to 

evaluate  the  integral  y 

p  +  x  =  4  r     du  (6i) 

c     J  o     V    u  +  ug(u) 


where 

a 

g    =  TEC 


o 


(^  -«o- 


When  the  dielectric  constant  varies  in  an  arbitrary,   but  montonic,   way 
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with  height,   the  integral  can  be  evaluated  readily  by  numerical 
methods.     Some  simplicity  is  obtained  if  it  can  be  assumed  that  the 
variation  is  exponential  in  form.     For  example  if 

K  -  Kq    =    (Kj/c)    (1  -  e"ch)  (62) 

it  is  seen  that  K  decays  exponentially  from  the  value  K     at  the  surface 
(h=0).     The  constant  c  is  of  the  order  0.13  for  an  average  atmosphere 

and  c  measures  of  the  rate  of  decay  of  the  dielectric  constant  with 
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increasing  height.     Before  proceeding  further  in  this  direction,   it  is 

probably  preferable  to  show  the  connection  of  the  above  results  with 

those  of  geometrical  optics. 

5.     Distance  to  the  Horizon 

The  distance  s  to  the  horizon  from  an  elevated  point  in  the 
atmosphere  is  often  used  conveniently  as  a  parameter  in  comparing 
one  type  of  vertical  refractive  index  profile  with  another.     Using  a  method 
based      on  a  stationary  phase  principle.it  readily  follows  from  Bremmer^ 
(and  others)  that 

h 

b  =  n0  c  -»l       a  (63) 

Jo      r  ^[rn(r)]»  -  n.a* 

where  r  =  a  +  h,   n(r)  is  the  refractive  which  is  a  function  of  r  and  n0 
is  the  surface  value  [  i.  e.   ntt  =  n(a)]  . 


In  the  case  of  exponential  atmosphere 
so  that 


n  =  1  +  be"ch    =    4K  (64) 


n.  =n]h=Q    =    1        ,    =    sfKQ 

The  term  under  the  radical  cai.  oe  written 
a 


a 


^a2 


^a* 


(■  •  ■"-")■  (■  *  ¥)'  -  (>  * ')  ] 

] 


1  +  2be"ch+    2-5-    -    1     -    2b 
a 


i*.    +    2b    fe'ch    -     1 


under  the  assumption  that  b  <  <  1  and  h/a  <  <  1.     To  the  same 
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n                                                                                           1 
approximation    — -    in  the  integrand  can  be  replaced  by    .     Consequently, 

the  integral  becomes 


11  o 

-     J  -  g  (65) 

°    \/(2h/a)    +    2b  (e"ch    -    1) 


which  is  equivalent  to  equation  (61)  derived  from  the  full  theory  (with  p  =  0). 
Now,    it  is  convenient  to  break  the  integration  into  three  parts,    such  that 

h  h»  h"  h   =h'" 


(66) 


o  ^h»  h 

and  thus 

s    =    s'     +    s"    +    sIM 
In  the  first  range  (0  =  h  =  s') 

h» 


£ 


dh 


(67) 


^J    (2h/a)    +    2b  Ach)  +  Sr-^ 


3 

which  is  valid  if  h1  is  chosen  such  that  the  (ch)    /6  term  in  the  expansion 

-ch 
of  e  is  negligible  compared  to  ch.      This  is  satisfied  if 

(ch1)  2/6    <  <  1 

(For  example  if,    ch1  =    -j-     then    (ch1)  2/6    =    l/24)     . 

Writing  h' 

...  f  dh 


N/(2h/a)-  21   '     >0        n/KTI  +  ah) 


(68) 


we  find 


where 


G      = 


\T2/a  -  2bc       n/g 

be2 
(2/a)     -     2bc 


1     /  ch' 

\l        1    +    Gh« 


-F-        tanil       s/      1    x   kl  (69) 
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If  the  effective  earth's  radius  is  defined  by 

1         _     1       .       K« 


where 


e 


K      =    (1  +  b)2    &    1  +  2b 


and  ,      y  i  v  2 

K-     =     ,*    6+b." 


h  =  0 


and  thus 


s    -2bc 


J_    =    J_    .    be 
a  a 

e 


so  that 

1 

\[Za/a 


«'  =  tt4^^  ^"J-rnip  <70> 


where 

c  /J_  J_ 

N^7a"e    ^a       "   ae 

if    ah'    <  <     1 

s«     =*  *  -s/  ah'  =    n/  2a  h'  (71) 

V  2 a/a  e 

'     e 

which  is  the  usual  approximation  for  the  distance  to  the  horizon  from 
a  height  h'  in  a  linear  atmosphere.     The  hyperbolic  tangent  formula 
above  is  actually  the  distance  to  the  horizon  from  a  height  h1  in  an 
atmosphere  defined  by 

n    =    1    +    b  (1  -  ch    +    (ch)2/2)  (72) 

for    0    =    h    =    h'.     This  is  an  approximation  to 

.      ,     ,     -ch 
n    =    1    +    be 

as  mentioned  above. 

We  shall  turn  our  attention  to  the  third  portion    (i.e.  h"5h=h'"=h0); 


'h»«    \/(2h/a)    +    2b  (e"ch  ^1) 
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Now  the  term  in  the  radical  is  written  in  the  following  form 

-ch 
(2h/a    -     2b)    +    2be-ch    =    (2h/a    -    2b)  ["l    +     ^/a  -  2b  1    '       (74) 

Letting  u  =  (2h/a  -  2b)  it  follows  that 


1 


=     J    An  (2b»n 


kl 


,n 


A  (2bT     e-chn 

-chN«-i  I           j —       A                         n  n 

2bc          \      2            ^u           n    i    ■>    *  u 

uf  1  +   )  n=0,  1,  2,  3.  .  . 


-cnba 
e  x 


n 


cna 


u 


u 


n+1/2 


(75) 


where  A       =1,     A,     = 
o  1 


1         .       _      1.3 
T'     A2    ~   ~2Ti" 


a  1-3-5  A  1-3-5-7  , 

A-    =    -    -* — a — 7-    *     A.     =    t — ^ — r — a —       ana  so  on 
3  2«  4»  6  4  2 .  4 •  6 •  8 


provided 


2b  e 


-ch 


u 


<    1  .     This  is  satisfied  if  2b  <  h/a    ^    h"/a    or  if 


-4 


h"  >  2ba     s     2x3x10       x  6360    ^   3.  9  km  for  any  value  of  c. 

-ch 
Actually,   this  is  more  stringent  than  necessary  since  e  <  1  . 

„,  cna 

Thus  u0      -  — £—  u 

s'"     =   Y  AJ2b)ne-cnba     C         e_T/,      du(a/2). 


n 


n=0, 1,2... 


u"      u 


E+T7T 


Now  consider  the  integral  which  occurs  in  the  above  expression, 

00 


Vm(pu)=    j" 


u         u 


:-pu 
m+1/2 


du 


(76) 


(77) 


and  integrate  it  by  parts    successively  to  yield  the  recurrence  relations 


Vm(P«)     =    (m  -  1/2) 
Vm-l<Pu»  "    (5  -  3/2) 


-pu 


m 


u 


e   r 


n/r 


-  p  v 


m 


m 


ZJ/T 


-    P    V„    ,  (P 


L  U 

and  so  on,    until 


lu)l 


(78) 
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where 


where 


V2(PU»    =   TV2)[^     ■    PV1(PU) 

00  _  _  1 

Vl  ((3u)   =    I      ^-TfT     du   =  (-J-)*  erfc  <^> 


(79) 


oo 


erfc  (n/^u)    = 


u        u 

2      I    ^ 

^/pu 


nTtF 


dx 


(80) 


is  the  complement  of  the  error  integral  which  is  extensively  tabulated." 
Finally 

(81) 
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i»i   = 


T 


^An(2b)ne-nbca  ["V^pu")  -  Vn(puJ 
n 


where 


(3    =     SJ^l   ,     U"     =    2  (—   -  b) 


and 


u0     =    2(^2    -b) 


It  would  appear  that  in  most  cases  h*  can  be  taken  equal  to  h"  and 
chosen  so  that  h1  =  h"  =  4  km  which  is  a  convenient  upper  limit  for  the 
quadratic  profile  and  a  lower  limit  for  the  expansion  formula.     Thus 

s  =  s'  +  s,n. 

—ch 
An  alternative  is  to  return  to  the  integral  s"  and  expand  e  to 

terms  of  third  order.     This  leads  to 

h" 
1  dh 


:ii    = 


where 


N/(2/a)  -  2bc 
be2 


Y  %7I 


N/h  (1  +  Yh)(l  +    6h) 
be 


(82) 


(F^y 


and      6    = 


if-r  "  bc 


This  integral  is  expressible  in  terms  of  elliptic  integrals  (see  Pierce 
No. 53 5  for  example). 


14 
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6.     Concluding  Remarks 

The  present  analysis  provides  a  strong  justification  for  the  validity 
of  an  equivalent  earth  radius  concept  when  the  refractive  index  of  the 
atmosphere  varies  in  a  linear  manner  with  height  in  agreement  with 
the  conclusions  of  Bremmer    ,   Norton      ,    Miller     and  Northover 
The  effect  of  the  non-linearity  becomes  important  at  great  heights  as 
pointed  out  by  Millington        and  others       '         who  have  suggested  that 
the  theory  for  a  linear  atmosphere  could  be  modified  by  simply  re- 
locating the  optical  horizon  point.     In  the  present  paper  a  complete 
mathematical  justification  is  given  for  this  step  provided  that  the 
variation  of  the  refractive  index  with  height  is  monotonic. 

The  idea  that  the  structure  of  the  diffraction  fields  in  a  monotonically 

varying  or  tapered  atmosphere  is  similar  to  that  for  a  homogeneous  or 
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airless  atmosphere  is  in  apparent  contradiction  to  Carroll's  theory 

of  "radio  twilight".     However,    the  model  he  uses  is  not  smoothly 
tapered  since  there  are  discontinuities  in  the  gradients  of  the  refractive 
index.     In  the  bilinear  model,    for  example,    the  atmosphere  is  represented 
by  a  layer  of  dielectric  with  a  constant  gradient  up  to  about  30,000  feet 

after  which  the  gradient  is  zero.    Carroll  obtains  his  solution  by  matching 

19 
the  appropriate  wave  functions  at  the  atmosphere-free  space  interface 

In  the  free  space  region  (above  30,000  feet  for  the  bilinear  model)  only 
outgoing  waves  are  allowed  but  at  the  top  of  the  dielectric  graded  layer 
the  two  types  of  wave  functions  which  are  used  can  be  identified  as  up- 
going  and  downcoming  waves.     In  the  present  treatment  which  uses  a 
modified  W.K.B.    method,   the  solution  in  this  region  does  not  have  this 
double  character. 

Other  models  used  by  Carroll  and  Ring       '         such  as  the  inverse- 
square  and  trilinear  profiles  can  be  objected  to  on  similar  grounds. 
It  is  of  interest  to  note  that  the  inverse-square  profile  is   somewhat 
better  in  that  the  discontinuity  is  in  the  second  derivative  at  the  bounding 
surface  and  the  computed  attenuation,    as  expected,    is  greater.      The 
trilinear  model  has  two  piece-wise  linear   sections  with  interfaces  at 
15,000  and  90,000  feet  and  the  mode  computations  show  less  attenuation 
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which  is  probably  a  consequence  of  the  strong  reflections  from  the 
lower  interface. 

It  is  the  opinion  of  this  writer  that  the  Carroll  and  Ring  computations 
are  correct  for  the  models  specified.     However,   they  cannot  be  regarded 
as  valid  solutions  for  propagation  in  a  monotonically  tapered  layer.     It 
is  probable  if  the  model  were  chosen  so  that  the  higher-order  derivatives 
of  the  refractive  index  are  continuous  at  the  interfaces,  the  mode  solutions 
and  the  modified  W.K.B.    solutions  should  approach  one  another.     Un- 
fortunately,  the  poor  convergence  of  the  resulting  mode  solutions  will 
make  this  ultimate  comparison  very  difficult. 
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Radio  Propagation  Engineering.  Data  Reduction  Instrumentation.  Modulation  Research.  Radio  Noise. 
Tropospheric  Measurements.  Tropospheric  Analysis.  Propagation  Obstacles  Engineering.  Radio-Meteor- 
ology.    Lower  Atmosphere  Physics. 

Radio  Standards.  High  Frequency  Electrical  Standards.  Radio  Broadcast  Service.  High  Frequency  Im- 
pedance Standards.     Electronic  Calibration  Center.     Microwave  Physics.     Microwave  Circuit  Standards. 

Radio  Communication  and  Systems.  Low  Frequency  and  Very  Low  Frequency  Research.  High  Fre- 
quency and  Very  High  Frequency  Research.  Ultra  High  Frequency  and  Super  High  Frequency  Research.  Modula- 
tion Research.     Antenna  Research.     Navigation  Systems.     Systems  Analysis.     Field  Operations. 


